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Abstract 



Based on density functional theory the influence of curvature on the de- 
pletion potential of a single big hard sphere immersed in a fluid of small hard 
spheres with packing fraction ij s either inside or outside of a hard spherical 
cavity of radius R c is calculated. The relevant features of this potential are 
analyzed as function of r] s and R c . There is a very slow convergence towards 
the flat wall limit R c — > oo. Our results allow us to discuss the strength of 
depletion forces acting near membranes both in normal and lateral directions 
and to make contact with recent experimental results. 
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Many biological processes are controlled by the interactions of macromolecules with 
cell membranes. Besides highly specific interactions of steric and chemical nature there 
are also entropic force fields which are omnipresent but whose actions depend only on gross 
geometrical features. These so-called depletion forces arise because both the membranes and 
the macromolecules generate excluded volumes for the small particles forming the solvent. 
Although these forces have been discussed for biological systems for many years IJ, the 
simultaneous presence of many other forces severely impedes the precise analysis of depletion 
forces in such systems. Therefore it is highly welcome that dissolved colloidal particles 
can be tailored such that they resemble closely the effective model of monodispersed hard 
spheres confined by hard walls 0. This allows one to study the depletion forces exclusively 
and to compare them quantitatively with theoretical results. Once a satisfactory level of 
understanding has been reached for these model systems one can apply with confidence this 
knowledge to the much more complex biological systems. Moreover, in colloidal suspensions 
themselves depletion forces can be exploited to organize self-assembled structures [[|. This 
dedicated use requires a detailed knowledge of the mechanisms of depletion forces, too. 

Along these lines in recent years there has been significant progress in understanding the 
depletion forces between two big spheres and between a single big sphere and a flat wall 
based on experiments [ffl-0], analytic results |I0-O|, and simulations [H. A generic feature 
of membranes is, however, that they are not flat. The variation of the local curvature leads 
to a new quality of the depletion forces in that they are no longer directed only normal 
to the surface, as for a flat wall or a wall with constant curvature, but that there is also 
a lateral component which promotes transport along the membrane. So far there are no 
systematic theoretical studies available which accurately predict this important curvature 
dependence of the depletion forces. Based on the quantitatively reliable density functional 
theory developed by Rosenfeld |T5|] we compute the depletion potential for a big sphere of 
radius Rb either inside or outside of a spherical cavity with radius R c , filled with a solvent of 
hard spheres of radius R s , as function of R c and of the packing fraction of the small spheres 
77 s defined as the fraction of the total volume occupied by the small spheres |16j . This enables 
us to make contact with a recent experiment in which the curvature dependence of depletion 
forces has been investigated by monitoring colloidal particles enclosed in vesicles [|I7J . 

Figure [l] shows our results for the depletion potential W(r; R c ) in units of k B T = [3~ l 
inside and outside of a hard spherical cavity centered at r = 0. Here W(r; R c ) is the 
difference of the grand canonical free energies of the hard sphere solvent in the presence and 
absence, respectively, of a big hard sphere whose center is fixed at a distance r from the 
center of the cavity. The choice of the parameters corresponds to those of the experiment in 
Ref. P7| . The solid line denotes the actual depletion potential compared with its so-called 
Asakura-Oosawa approximation (AOa) |18|] (dotted line) which is valid only for such small 
values of r] s that the solvent can be treated as an ideal gas and which is determined by 
the overlap of the excluded volumes around the big sphere and the hard wall. The AOa, 
which has been used to interpret the experimental data in Ref. Hl7| , satisfactorily predicts 
the value W c of the depletion potential at contact and its derivation with respect to R c , but 
otherwise it fails considerably. Whereas WAOaij] R c ) is purely attractive and vanishes for 
r R c ^f Rb^f2R s , the actual potential is both attractive and repulsive with the wavelength 
of the oscillations approximately given by 2R S . The correlation effects of the solvent generate 
the potential barrier AW r (see Fig. [I]) which is completely missing within the AOa, and they 
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increase the range of the potential significantly beyond the AOa range of 2R S . The presence 
of this potential barrier has very pronounced repercussions for the diffusion dynamics of the 
big particle. The time r required to overcome a potential barrier AE is proportional to 
exp((3AE). Therefore, with AE = AW r , for the big sphere it takes about e 2 ' 5 ~ 12 times 
longer to reach the wall from the center of the cavity or to escape from the wall as compared 
with the time estimated from the AOa (see Fig. p]). This difference still awaits experimental 
confirmation. For larger packing fractions r] s this difference becomes even larger. Thus the 
value W c = AW r — AW e of the depletion potential at contact can be obtained from the 
ratio of rates of escaping from the wall and reaching it, respectively. Taking into account 
only the former one, as suggested by the AOa, would be misleading. The overall structure 
of the depletion potential outside of the cavity is similar to the one inside, but its amplitude 
is reduced (Fig. [1]). This is in line with the trend predicted by the AOa which expresses 
the fact that the overlap between the excluded volumes around a big sphere and the wall 
is larger inside the spherical cavity than outside. This geometrical difference between the 
excluded overlap volume inside and outside of the cavity is schematically drawn as inset in 
Fig. [I]. In the inset the boundaries of the excluded volumes around the big sphere and the 
cavity wall are indicated by dashed lines and the overlap of excluded volumes is shaded. 

In Figs. |2] and |^ we discuss as function of the packing fraction rj s the relevance of the cavity 
curvature for the depletion potential relative to the case of a planar wall, i.e., R c — ► oo by 
focusing on the potential at contact W c and the escape potential barrier AW e , respectively. 
As expected in the limit R c — > oo both quantities reach a common value for the outside 
and inside potential. These common values are in very good agreement with simulation 
data for the planar wall [I|J]. For all values of R c the absolute values of the potential 
parameters inside are lager than outside. They increase stronger than linearly with increasing 
i] s . For increasing curvature 1/R C the outside potential becomes weaker whereas the inside 
potential becomes much stronger; this holds for all values of rj s . The gain of strength of 
the inside potential upon increasing the cavity curvature is much more pronounced than the 
corresponding loss of strength on the outside. This difference between the behavior inside 
and outside widens strongly with increasing rj s . The dependence of the depletion potential 
on the cavity curvature is surprisingly strong even for large values of R c /R s . For rj s = 0.3 
and R c = 50R S the density profile p s (r) of the solvent without the big sphere near the curved 
wall differs only slightly from that near a flat wall; its contact value differs by 1.5% from 
the corresponding flat one. However, the potential at contact still differs by almost 10% 
from the corresponding flat value. This amplification of the influence of the curvature can 
be understood in terms of the geometric considerations of the overlap volumes leading to 

pwrm - -n, (i + ^) (i) 

with s = Rb/R s and 7 = +1 and —1 for outside and inside, respectively. Within AOa 
the absolute value of the amplitude of the first curvature correction is the same inside and 
outside as can be seen from 

p W AOa {oQ) " 1 " 3S + 1 RJR S + W Rs) )• (2) 

For R b ^> R s , i.e., s ^> 1 the curvature dependence of W c is significantly enhanced. 
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In a recent experiment |17j Dinsmore et al. used video microscopy to monitor the position 



of a single big colloid particle with = 5.7 R s immersed in a solution of small colloids 
with r] s = 0.3 inside a vesicle. These were charge stabilized colloids which are supposed 
1% to resemble closely the model system of hard spheres confined by a hard wall. From 



the probability p(x) of finding the big colloid at the position x - the quantity actually 
measured in the experiment - one can infer the depletion potential W(x) according to p(x) = 
Pbuik^~^^ with W(x) = in the bulk. Thus the spatial resolution of the experimentally 
determined W(x) reflects that of p(x), which at best was approximately 2R S in Ref. [|T7| . 



In the experiment the big colloid particle was observed to be most of the time very close to 
the vesicle wall inside a shell whose width was chosen to be 6.7R S , i.e., even larger than the 
optimal resolution 2R S . Approximating W(x) by W(r; R c ), where r is the minimal distance 
between x and the vesicle wall and R c the local radius of the wall at this closest point 
(see below), the experimentally determined depletion potential therefore corresponds to an 
average of W(r; R c ) shown in Fig. [I] over all visible oscillations. Within this shell the big 
colloid was observed to be more often in regions with a small local radius of curvature. This 
is in qualitative agreement with both the AOa and our present approach. Quantitatively, in 
Ref. |17| the ratio p s heii/Pbuik of the probabilities of finding the big particle within the shell 



and in the bulk, respectively, was determined. The authors have documented the logarithm 
of this ration, denoted as 

(3AF := - In ^ ~ - In ( f d 3 re^ w ^ /V shel ) , (3) 

Pbulk \Jv ahell J 

relative to its value flAF^ near a flat wall, i.e., for R c — > oo; V s h e ii is the volume of the 
shell. In Fig. ^ we compare the published data for f3AF — flAF^ with our prediction for this 
quantity based on the actual depletion potential as shown in Fig. [I], with the corresponding 
prediction if in Eq. @ and for /SAF^ the full AOa is used, and with the so-called truncated 
AOa, f3W^° a -(3W^° a , i.e., the AOa values for W(r; R c ) at wall contact, which was used for 



comparison with theory in Ref. [[L7]]. We find that our theoretical prediction as well as the full 
AOa are closer to the experimental data than the truncated AOa. Given the high accuracy of 
our calculation the remaining discrepancy cannot be due to insufficient theoretical knowledge 



of W(r;R c ). The fact that in Ref. |T7| only the radius of curvature in plane could be 
determined, polydispersity of the solvent, and the possible presence of dispersion forces are 
among the candidates for explaining this discrepancy. The small difference between the full 
AOa and our theoretical results in Fig. |] demonstrates that the present spatial resolution 
cannot discriminate between the rich structure shown in Fig. [I] and its AOa. This emphasizes 
the need for future experiments with significantly increased spatial resolution. 

For a cavity of general shape z = /(R = (R x ,R y )) of its surface relative to a suitable 
chosen (x, y) reference plane our results allow us to determine approximately lateral depletion 
forces. To this end we introduce normal coordinates |19| (R x , R y , r) such that x = (x, y, z) = 
(R x , R y , /(R)) + rn(R) where r is the minimal distance of the point x from the surface /(R), 
(R x ,R y ) (7^ {x,y)) are the lateral coordinates of that point on the surface closest to x, and 
n = (— V/(R), 1)/ a/1 + (V/) 2 is the local surface normal pointing towards x. The actual 
depletion potential W(x; [/]), which depends functionally on /(R), can be approximated 
by W(x; [/]) ~ W(r(x); i? c (R(x))). Here we have assumed that the cavity surface varies 
sufficiently smoothly so that its local two principal radii of curvature can be described by 
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the single radius R C (H). Within this approximation the components iy a 4 t (x) of the lateral 
depletion force in the direction of the tangential vector tj, i = 1, 2, is given by 



with r = r(x; [/]), R = R(x; [/]), V = d/dx, and n • t 4 = = ti • t 2 . The factor dW/dR c 
can be determined from our above results. Its absolute value increases with decreasing 
radius of curvature. Its sign differs inside and outside of the cavity and changes as func- 
tion of the normal distance r. A big sphere approaching the convex nonspherical wall from 
inside (outside) is exposed to an oscillating lateral force, which pulls the sphere towards 
regions of the surface with a small local radius of curvature close to a minimum (maxi- 
mum) of the depletion potential and pushes it away from these regions close to a maximum 
(minimum). These oscillations of the lateral force originate from the packing effects of the 
small spheres. When the big sphere has reached the wall the lateral force pulls it along 
the surface to that point with the smallest local radius of curvature inside the cavity and 
pushes it away from this point towards the point with the largest local radius of curvature 
outside of the cavity, respectively. This can already be understood within the AOa which 
predicts d((3W^° a (R c ))/dR c = -3rj s -fsR b (R c + ^R h )- 2 (see Eq. (§). Since Eq. flip provides 
a very good description of the contact value of the depletion potential, even for high packing 
fractions, the above expression for its derivative is also quantitatively reliable. 

Applying Eq. (f|) for an ellipsoidal cavity with semi-axes a = b = 20R S and c = 30R S and 
for a big sphere with R b = 5R S immersed in a fluid of small spheres with rj s = 0.3 yields at 
contact a maximal ratio of the lateral force to the normal force of \Fi at / F norm \ = 0.07 inside 
and outside of the cavity. Such quantitative predictions for the size of entropic lateral forces 
along curved confining walls still await tests by simulations and experiments. 
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FIGURES 




FIG. 1. Depletion potential W(r; R c ) in units of ksT = (3~ l inside and outside of a spherical 
cavity centered at r = with radius R c = 4QR S (full vertical line). The radius of the big sphere is 
Rb = 5.7R S and the packing fraction of the hard sphere solvent (radius R s ) is r] s = 0.3. Due to the 
hard core interactions the center of the big sphere cannot enter the space indicated by the dashed 
lines. The solid line is obtained from density functional theory whereas the dotted line corresponds 
to the Asakura-Oosawa approximation which is valid for small values of i] s and exact in the limit 
i] s — > 0. W c < 0, AW r > 0, and AW e = —W c + AW r > denote the depletion potential at contact 
with the wall and the potential barriers for reaching the wall and escaping from it, respectively. 
The inset shows the different geometry of the (shaded) overlap of excluded volumes of a big sphere 
(hatched) inside and outside of a spherical cavity (hatched). (For reasons of clarity the schematic 
drawings correspond to Rb = 0.7R S and R c = 8.6R S .) 
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FIG. 2. Depletion Potential at contact 0W C (see Fig. ||) as function of the cavity curvature 
(R c /R s )~ l for three packing fractions rj s outside (upper curves) and inside (lower curves) of the 
cavity with = 5R S . The symbols at R~ l = indicate simulation data for a flat wall Jl4[ . 
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FIG. 3. Escape potential barrier (3AW e for a big sphere with radius Rb = 5R S inside (upper 
curves) and outside (lower curves) of a spherical cavity of radius R c for three values of r] s . The 
symbols at R~ l = indicate simulation data for a flat wall [14]. 
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FIG. 4. Comparison of the experimental data [17] for (3AF — PAF^ (Eq. (|3|)) as function of 
the local radius of curvature R c in units of R s with the prediction of the present theory (full line in 
Fig. |]) for this quantity. The dotted line corresponds to the dotted line in Fig. |l]. The dashed line 
denotes f3W^° a — (3W^^ a which are the AOa values for W(r; R c ) at wall contact which served in 
Ref. |17| as theoretical comparison. 
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